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Tables for Use in the Interpretation of Paramagnetic 
Behavior Below 1° K; for the Chromic Alums (J = 3/2) 

Ralph P. Hudson and Charles K. McLane 

The method of interpreting susceptibility and entropy data, obtained in adiabatic 
demagnetization experiments, by means of the theory of Hebb and Purcell in order to derive 
absolute temperatures is recounted briefly. This forms an introduction to tables that 
enable one to carry out this procedure for the case of the chromic alums. 



The low- temperature behavior of the chromic 
alums is largely determined by the fact that the 
fourfold degenerate ground level of the Cr +++ ion is 
split by the crystalline electric field (Stark effect) 
into two doublets separated by an energy 5. The 
theory of Hebb and Purcell [1] * yields expressions 
for the susceptibility, xo, and the entropy, S, as 
functions of the absolute temperature, T, and this 
splitting 8. Thus one may obtain a relation between 
S and xo for any chosen 5. In practice one usually 
employs the "magnetic temperature" T 7 * instead of 
Xo, and T*=cJxo, where c is the Curie constant per 
cubic centimeter. 

Now experiment also leads to a xo versus S (or 
T* versus S) relation, for S may be calculated quite 
accurately, simply from the conditions of magnetic 
field and temperature at the beginning of the isen- 
tropic demagnetization process, and xo is measured 
at the end. One therefore has only to find that value 
of 5 which leads to an S versus T* curve fitting the 
experimental points. Then over the range of agree- 
ment one may suppose the Ilebb and Purcell xo(T) 
and S(T) formulas to be individually correct, and 
the former leads at once to the establishment of 
the absolute temperatures. 

Table 1 gives, effectively, corresponding values of 
T/8 and T*/8 calculated from the Hebb and Purcell 



Table 1. Corresponding values of kT/8 and k(T*—T)/8 cal- 
culated from the Hebb and Purcell [1] formula: T*/8 = y~ 1 X T/8 
(see text) 





T*-T 




T *_ T 




T*-T 


kT/8 


8/k 


kT/8 


8/k 


kT/8 


8/k 


0.05 


0.0281s 


1.0 


0. 0312 5 


3.0 


0.01103 


.10 


.04705 


1.1 


. 0287 2 


3.2 


.0103s 


.15 


.05842 


1.2 


. 0265 5 


3.4 


. 00974 


.20 


. 0639 


1.3 


.02467 


3.6 


.0092 2 


.25 


. 0653 4 


1.4 


.0230s 


3.8 


.0087i 


.30 


.0643i 


1.5 


.02159 


4.0 


.0083i 


.35 


.0619? 


1.6 


. 0203o 


4.2 


.0079o 


.40 


.0590i 


1.7 


.01917 


4.4 


. 0075 4 


.45 


.0558s 


1.8 


.01814 


4.6 


.00722 


.50 


.0527l 


1.9 


.01722 


4.8 


.0069 3 


.55 


.0497i 


2.0 


•0163 9 


5.0 


.OO664 


.60 


.0469 L 


2.1 


.0156 3 






.65 


. 0443 2 


2.2 


. 0149 4 






.70 


.0419$ 


2.3 


.0143i 






.75 


. 0397 5 


2.4 


.0137s 






.80 


. 0377 5 


2.5 


.0131s 






.85 


. 0359 L 


2.6 


.0127o 






.90 


. 0342s 


2.7 


.0122s 






.95 


. 0326s 


2.8 


.0118o 






1.00 


. 0312 5 


2.9 
3.0 


.01139 

.01103 







formula [1]: 7 7 */5=^ 1 X7 7 /5, where 



1 Figures in brackets indicate the literature references at the end of this paper. 
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Tables 2 to 5 contain the data necessary for cal- 
culating the "experimental entropy/' that is, in the 
magnetic field at the initial temperature, and the 
field-free entropy as a function of kT/5. 

Considering first the case of the entropy in zero 
field and at low temperatures, two effects contribute 
to a reduction in the entropy below the "ideal" 
value of R log e 4, namely, the Stark splitting and 
magnetic interaction. The theory treats these two 
contributions separately, and they are additive. 
The "Stark entropy" is calculated from the partition 
function, Z c = 2(l-{-e~ d/kT ), and is given in table 2, 
Table 3 lists as a function of x (where x=kT/8), 
the quantity (d/dx)(Q/x)j 12 being as defined in 



reference [1]. 
given by 



The magnetic entropy, S m , is then 



s m /B=(r/sy 



d_ 
dx 



(!) 



(2) 



where t/* = 3c(=0.0211 deg and 0.0194 deg for the 
potassium and methylammonium alums, respec- 
tively) . Hence 



log e 4- 



-S/R= (log, 4-S/£) stark + SJIi. 



Table 2. "Stark entropy", or the reduction in entropy below 
the value R log c 4, caused by the effect of the crystal field 



kT/8 


log ei-S/R 


kT/8 


log, 4- S/R 


kT/8 


\og e 4-S/R 





0. 6931 5 


1.0 


0. 1109 4 


3. 


0. 0137o 


0.05 


. 6931 5 


1.1 


. 0935a 


3.2 


. 0120e 


.10 


. 6926 5 


1.2 


.0798i 


3.4 


. 0107 


.15 


. 6834 


1.3 


.O6883 


3.6 


. 0095s 


.20 


. 6529 7 


1.4 


. 0599 3 


3.8 


. 0085 8 


.25 


. 6030s 


1.5 


.05262 


4.0 


. 0077s 


.30 


. 5432 8 


1.6 


.04654 


4.2 


. 00704 


.35 


. 4821 2 


1.7 


. 0414 5 


4.4 


. 0064 2 


.40 


.4246i 


1.8 


.03714 


4.6 


. 0058 7 


.45 


. 3729 9 


1.9 


. 0334 6 


4.8 


.0054o 


.50 


.3278i 


2.0 


. 0303o 


5.0 


. 0049 8 


.55 


.2888i 


2.1 


. 0275 6 






.60 


. 2553 6 


2.2 


. 0251 7 






.65 


. 2267o 


2.3 


. 0230 8 






.70 


. 20214 


2.4 


. 02124 






.75 


. I8IO4 


2.5 


. 0196i 






.80 


. 1628 4 


2.6 


,0181s 






.85 


. 14709 


2.7 


.01686 






.90 


. 1333 9 


2.8 


.0156 9 






.95 


. 1214 3 


2.9 


.0146s 






1.00 


. 1109 4 


3.0 


.0137o 
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Table 3. Corresponding valves of kT/d and the function 
(d/dx)(n/x) [1] 

The "magnetic entropy", S m , or the entropy diminution additional to the 
Stark entropy (table 2), is obtained from the relation S m /R= (j/8) i (d/dx){Q,lx) , 
(see text) . 



kT/8 


d /a\ 
dx \x J 


kT/8 


d /n\ 
dx \x J 


kT/8 


dx \x J 







1.0 


0. 7706 6 


3.0 


0. 1228s 


0.05 
.10 




1.1 
1.2 


. 6795i 
. 6005 6 


3.2 
3.4 


. 1090 2 
.09715 


7. 1616 6 


.15 


2. 9618 7 


1.3 


. 5327i 


3.6 


. 0871o 


.20 


1.4571 e 


1.4 


. 4745 3 


3.8 


. 0785 2 


.25 


0. 9478 8 


1.5 


. 4245 9 


4.0 


.07114 


.30 


. 8777 7 


1.6 


.38i7o 


4.2 


. 06474 


.35 


.9611s 


1.7 


. 3444 5 


4.4 


.05917 


.40 


1. 0705 6 


1.8 


.3122i 


4.6 


. 0542 9 


.45 


1. 1518 7 


1.9 


. 2841 1 


4.8 


. 0499 8 


.50 


1. 1940e 


2.0 


. 2595o 


5.0 


. 0461 6 


.55 


1. 2007 6 


2.1 


. 2378 5 






.60 


1. 1820e 


2.2 


. 2187 6 






.65 


1. 1451 3 


2.3 


. 2017 5 






.70 


1.09676 


2.4 


. I8667 






.75 


1.04258 


2.5 


. 17313 






.80 


0.98575 


2.6 


. I6IO0 






.85 


.9292o 


2.7 


. I0OO9 






.90 


. 8738 6 


2.8 


. 1402 3 






.95 


. 8208 5 


2.9 


. 1313o 






1.00 


. 7706 6 


3.0 


. 12283 







Table 4 lists values of the entropy, for correspond- 
ing values of H/T, calculated by means of the 
"Brillouin formula. " For the case «7=3/2, one has 



(S/R) B =acoiha — 4<z coth 4a +log e sinh 4& - 



-log e sinhu\ 
(3) 



where a—fiH/kT; fi is the Bohr magneton, 
9.271 X10 -21 emu, and k is the Boltzmann constant, 
1.3805X10" 16 erg deg" 1 . 

In order to derive an accurate value of 8, however, 
one must take cognizance of the fact that the 
Brillouin entropy (table 4) is only appropriate to the 
ideal case of a fourfold degenerate ground level. 
The effect of the crystalline electric field is to modify 
the pattern of the level splitting in a magnetic field, 
and for small fields the correction, AS/R, to be 



applied to the entropy reduction, log e 4 — S B /R, is 
relatively large. 

An entropy formula has been derived [3] that takes 
this effect into account, but it is more unwieldy than 
the Brillouin formula (a straightforward function of 
HjT) as it is a function of H, T, and 8. For the 
chromic alums the value of S/k is always close to 0.25 
deg, and starting temperatures (that is, bath temper- 
atures) will, in general, be of the order of 1.1° to 
1.2° K. Because AS/R<^S B /R, the Brillouin and 
exact entropy curves run very close together, and 
the following approximate method of calculating a 
new AS/R under slightly changed conditions of 8 
and/or T (say §1, I\), given a table of exact values 
for conditions <5 , T Q , is found to be very useful. 

Let A=5/2, a=A/T, and = /x#/A. One may then 
take the expressions for both the Brillouin entropy 
and the exact entropy and obtain their derivatives 
with respect to a. Denoting derivatives by a 
primed symbol, we then have 

(S B IR) a ~(S B /R) ao +(S' B /R) ao X(a-a ), (Brillouin) (4) 
and 
(S/R) a « (S/K) ao +(S'/R) aQ X(a- ao). (exact) (5) 

Subtracting (5) from (4) — -for a positive quantity — 
we find 

(AS/R) a = (AS/R) aQ + (AS'/R) ao X(*- a ) (6) 

to a satisfactory degree of accuracy. It must be 
emphasized that eq (4) and (5) are only approximate. 
Higher derivatives have been neglected, but these 
effectively disappear in the subtraction, which leads 
to eq (6). 

In table 5 the entropy values have been calculated 
for various values of H, and with 8/k = 0.27 deg, 



Table 4. Entropy as a function of a( = /jiR/kT) calculated from the " Brillouin formula" : S/R = a coth a — 4a coth 4a + log c sinh 

4a— loge sinh a 

Tables for J=l/2, 3/2, 5/2, 7/2 can be found in an article by Hull and Hull [2]; the greater part of table 1 has been taken from this source. 



a 


S/R 


1 
log, 4- S/R 





S/R 


loge 4- S/R 


a 


SIR 


log. 4- S/R 





1. 38629 





0.40 


1. 07179 


0. 3145 


1.50 


0. 20817 


1.1781 


0.01 


1. 38604 


0. 0003 


.45 


1.00961 


. 3767 


1.55 


. 19228 


1. 1940 


.02 


1. 38530 


.0010 


.50 


0. 94754 


.4387 


1.60 


. 17756 


1. 2087 


.03 


1. 38405 


.0022 


.55 


. 88658 


.4997 


1.65 


. 16393 


1. 2224 


.04 


1. 38231 


.0040 


.60 


. 82750 


. 5588 


1.70 


. 15131 


1. 2350 


.05 


1. 38007 


.0062 


.65 


. 77083 


. 6155 


1.75 


. 13963 


1. 2467 


.06 


1. 37735 


.0089 


.70 


. 71690 


.6694 


1.80 


. 12882 


1.2575 


.07 


1. 37415 


.0121 


.75 


. 66591 


.7204 


1.85 


. 11882 


1.2675 


.08 


1. 37047 


.0158 


.80 


. 61796 


.7683 


1.90 


. 10957 


1. 2767 


.09 


1. 36632 


.0200 


.85 


. 57303 


.8133 


1.95 


.10102 


1.2853 


.10 


1.36171 


.0246 


.90 


. 53102 


.8553 


2.00 


.09311 


1. 2932 


.12 


1.35115 


. 0351 


.95 


. 49185 


.8944 


2.1 


. 07904 


1. 3072 


.11 


1. 33887 


.0474 


1.00 


. 45543 


. 9309 


2.2 


. 06705 


1. 3192 


.16 


1. 32496 


.0613 


1.05 


. 42156 


.9647 


2.3 


. 05681 


1. 3295 


.18 


1.30952 


.0768 


1.10 


.39011 


.9962 


2.4 


. 04809 


1. 3382 


.20 


1. 29266 


.0936 


1.15 


. 36091 


1. 0254 


2.5 


. 04068 


1. 3456 


.22 


1. 27449 


. 1118 


1.20 


. 33382 


1.0525 


2.6 


. 03438 


1.3519 


.24 


1. 25514 


.1311 


1.25 


. 30871 


1.0776 


2.7 


. 02903 


1. 3573 


.26 


1. 23471 


.1516 


1.30 


. 28543 


1. 1009 


2.8 


. 02449 


1.3618 


.28 


1. 21334 


.1730 


1.35 


. 26386 


1. 1224 


2.9 


. 02064 


1. 3656 


.30 


1. 19115 


.1951 


1.40 


. 24386 


1. 1424 


3.0 


.01739 


1.3689 


.32 


1. 16824 


.2181 


1.45 


. 22534 


1. 1609 








.34 


1. 14475 


.2415 


1.50 


. 20817 


1. 1781 








.36 


1. 12077 


.2655 














.38 


1. 09642 


.2899 














.40 


1. 07179 


.3145 
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!r=l.l K (the Brillouin entropy itself is, of course, 
independent of 8, but the values of H are listed in 
terms of tiH/A, which appears in the analytical 
expressions for the exact entropy). We thus have 
A = 0.K*5 and a =0. 1227273. The derivatives are 
given in columns 3 and 5. 

Now, in any experiment, T is determined by 
experimental conditions, and 8 is usually known 
fairly accurately. (It is not illogical to derive an 
accurate value of <5 from entropy data, the calcula- 
tion of which involves 8; it is essentially the method 
of successive approximation.) These together fix 
the value of a, and a tabulation of AS/R is readily 
obtained by making use of table 5 and eq (6) . 

One then plots a curve of AS/R versus (S B /R)a 
(using eq (4)) to obtain the correction curve appro- 

Table 5. Brillovin entropy, entropy correction, and first 
derivatives with respect to a( = A/T) for the "standard value" 
ao of 0.1227273 ( = 0.135/1.1), for varying values of 0( = /J7/ 
A). A = 6/2 (see text) 



fi-l*H/A 


(S B IR)a 


-(S'fi//?)« 


(AS(R) a , 


(AS'IR)a 





1. 38629 





0. 00747 


0. 1209 


0.4 


1. 38029 


0. 097397 


. 00744 


.1197 


.8 


1. 36259 


. 380323 


.00734 


.1162 


1.2 


1. 33400 


.821912 


. 00716 


.1110 


1.9 


1. 26186 


1 . 866246 


. 00674 


.0983 


2.8 


1. 14042 


3. 337499 


. 00605 


. 0790 


3.7 


1. 00*51 


L 610410 


. 00529 


. 0600 


4.6 


0. 86918 


5. 478282 


,00454 


.0135 


6.0 


. 67952 


6.036330 


. 00350 


. 0238 


7.5 


. 51466 


5.909666 


. 00260 


. 0097 


9.0 


. 38737 


5.417031 


. 00191 


.0010 


10.5 


. 29057 


4. 787. r .3 J 


.00137 


-. 0038 


12.0 


. 21734 


4. 128350 


. 00093 


-. 0022 



priate to one's own experimental conditions. For 
any one demagnetization, S B /R is readily obtained 
from the H/T value and table 1, and AS/R is read 
off from the correction curve. 

To illustrate the efficiency of the method, table 6 
lists values of AS/R calculated rigorously and by 
the above approximation for the case a=0.H3043 
(5-0.26 deg, 2 T =1.15° K) and for the case a= 
0.133333 (5 = 0.28 deg, T=1.05° K). Inasmuch as 
one is only interested in an accuracy to the fourth 
decimal place, or two significant figures, it may be 
seen that the method is very useful. 

Table 6. Comparison of entropy-correction values calculated 
exactly and by the approximate method for the cases (a) 
a = 0.1 13043 (8/k = 0.26 deg, T=1.15°K), and (b) a = 
0.133333 (5[k = 0.28 deg, T=1.0S°K) 



fi=liHlA 


a = 0.133333 


a = 0.1 13043 


AS/R 

(exact) 


ASfR 
(approx.) 


AS/R 

(exact 1 


AS/R 
(approx.) 



3.7 

10.5 


0. OO881 

.005'.), 
.00132 


0. 00876 
. 0059 2 
.00133 


0. 0063 5 
.0047o 
.0013s 


0. 0063o 

.0017,1 
. 0014i 



[1] M. H. Hcbb and E. M. Purcell, J. Chem. Phys. 5, 338 

(1937). 
[2] J. II. Hull and R. A. Hull, J. Chem. Phys. 9, 465 (1941). 
[3] R. P. Hudson, Phys. Rev. 88, 570 (1952). 

Washington, September 17, 1953. 



2S2952— 54- 



35 



